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Limits on R-parity violation from cosmic ray antiprotons 



OO 
On 

3 



(N 
> 



^1" 
O 

a: 

Oh! 
43 



- 1—1 

x 



Edward A. Baltz* 
Department of Physics, University of California, 
601 Campbell Hall, Berkeley, CA 94720-3411 

Paolo Gondolo^ 
Center for Particle Astrophysics, University of California, 
301 he Conte Hall, Berkeley, CA 94720-7304 
and 

Max- Planck- Institut fur Physik, Werner-Heisenberg-Institut, * 
Fdhringer Ring 6, 80805 Miinchen, Germany. 



We constrain the hadronic R-parity violating couplings in extensions to the minimal supersym- 
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I. INTRODUCTION 

In supersymmetric models, gauge invariance allows in- 
teractions that violate lepton- or baryon-number conser- 
vation. To avoid them, a discrete symmetry called Im- 
parity is often imposed. The R-parity of a particle is 
defined as R = (— i) 2S, + L + 3B ; where S is its spin, and L 
and B are its lepton and baryon numbers. R is +1 for 
all standard model particles and —1 for all superparti- 
cles. Alternatively, it is an interesting exercise to allow 
R-parity violating couplings and constrain them by ex- 
perimental observations or cosmological considerations. 
With R-parity violation, the lightest superpartner decays 
into standard model particles. We will assume that the 
lightest superpartner is the neutralino. 

Antiprotons are a rare component of cosmic rays, most 
probably produced in spallation reactions in the interstel- 
lar medium. The measured antiproton flux limits any 
additional antiproton production in the galaxy. We find 
a powerful bound on the lifetime of the lightest super- 
partner from the observed flux of cosmic ray antiprotons. 
This allows us to bound the R-parity violating interac- 
tions that appear in extensions to the minimal supersym- 
metric standard model. 

We find the relic density of the lightest superpartner 
(LSP) using the standard cosmological model. We then 
allow these particles to decay via the hadronic R-parity 
violating operators. We employ simple models for the 
production rate and the propagation of the antiprotons. 
We can limit the R-parity violating couplings with ob- 
servations of the antiproton spectrum in cosmic rays. 

Our results strictly apply when the lifetime of the neu- 
tralino is longer than the age of the universe. However, 
there are other cosmological limits based on nucleosyn- 
thesis in the early universe, recent structure formation, 
diffuse photon and neutrino fluxes, and spectral distor- 
tions in the cosmic microwave background that cover a 
large fraction of neutralino lifetimes shorter than the age 
of the universe, and may well completely close the gap 
between this result and laboratory constraints. 

II. SUPERSYMMETRIC MODEL 

The supersymmetric model we consider is a minimal 
extension of the standard model with explicit R-parity 
breaking, defined by the superpotential 

w = /..,//, - h[; i rfjji: + h%u?QjH 2 

-nHiH 2 + \' ijk LiQjDl + X'! jk UtD1D%. (1) 

i,j,k are generation indices; Li and Qi are the SU(2)- 
doublet lepton and quark superfields; Ef, Uf and D\ are 
the singlet superfields; and H\ and H 2 are the doublet 
Higgs superfields. h l J , h l J and h l J are Yukawa matrices, 
and \i is a free parameter [|lj . 

The last two terms in the superpotential violate lepton- 
and baryon-number conservation explicitly. For three 



generations, there are 27 A' couplings and 9 A" couplings, 
since A" jfe is antisymmetric under exchange of the last two 
indices. 

Supersymmetry is broken softly by a scalar potential 
that in general contains many unmeasured parameters. 
For our purposes, a simple parameterization suffices. We 
assign a common mass scale mo to the sleptons and 
squarks at the electroweak scale, we relate the gaugino 
mass parameters Mi, M 2 , M3 by means of GUT relations, 
and we keep the trilincar soft terms A t and only for 
stop and sbottom. Electroweak symmetry is broken at 
tree level by assigning vacuum expectation values v\ and 
v 2 to the neutral scalar components of H\ and H 2 . The 
ratio tan/3 = v 2 /vi and the pseudoscalar Higgs boson 
mass ttia remain as free parameters. The total number 
of parameters is therefore seven. Further details on the 
class of models we consider are given by Bergstrom and 
Gondolo @. 

We use the scans of supersymmetric parameter space 
in Bergstrom, Edsjo, and Gondolo ||. They consist of 
~ 10 4 models that satisfy the experimental constraints 
reported by the Particle Data Group Q] , plus the bounds 
on the Higgs sector and chargino masses from LEP (we 
update to LEP2W B) and the constraints on the b — > sj 
rate from CLEO f|. The scans span values of and 
\M 2 \ up to 5 TeV, of m up to 3 TeV, of \A t \ and \A b \ 
up to 3too, of rriA up to 1 TeV, and of tan/3 between 1 
and 50. This range extends even beyond the parameter 
region usually considered as plausible. 

III. COMPOSITION OF THE GALACTIC HALO 

We do not assume that the dark matter is supersym- 
metric in nature. Given a supersymmetric model, the 
relic density of neutralinos is fixed and calculable in the 
standard cosmological model. We use the relic density 
computed by Edsjo and Gondolo 0, who include all neu- 
tralino self-annihilation diagrams, resonance and thresh- 
old effects, and chargino-neutralino coannihilations. In 
every model that is cosmologically allowed, the lightest 
neutralino is a cold relic from the Big Bang. Being a cold 
relic, we expect that it is attracted into the gravitational 
potentials of galaxies, and is at present a component of 
the cold dark halo. 

We use the simple estimate that the dark matter den- 
sity of the galactic halo at the position of the solar system 
is 0.3 GeV cm~ 3 . Neutralinos may make up the entire 
halo or only some fraction /. This fraction / > 17 x /f2DM, 
where £Ia is the cosmological density of species A mea- 
sured in units of the critical density. 

We now need an estimate of JIdm- A common prescrip- 
tion H is to take f^DM — ^x if ^x^ 2 — 0.025, which is the 
value indicated by dynamical mass measurements of the 
amount of dark matter in galactic halos. If Vl x h? < 0.025, 
this prescription sets Qdm^ 2 = 0.025. This procedure is 
not conservative enough for our purposes because there 
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may be more than one cold dark matter component in 
the galactic halo. Therefore, we use the constraint that 
Udm'' 2 < 1, based on the age of the universe ||, to give 
a lower bound on / = Q, x h 2 . In this way we obtain a con- 
servative limit on the R-parity violating couplings. This 
may be overly conservative, and indeed our final results 
would improve by approximately a factor of five if we 
used the less conservative fraction. 



IV. DECAY RATE 

We consider neutralino decay through both the lepton- 
numbcr violating term X' LQD C and the baryon-numbcr 
violating term X"U C D C D C in the superpotential. The 
coupling A^ fe gives rise to the decays x ~ > v idjdk, X ~ * 
e^Ujdk and their charge conjugates. The coupling A^ fe 
gives rise to the decays x ~ * u idjdk (j ^ k) and their 
charge conjugates. 

The differential decay rates for these processes are cal- 
culated in Ref. ||. Here we quote the decay rate for 
X — > uds in the limit where the squarks are heavy and 
degenerate, 



T(x -> uds) 



32tt 



3 m x |Au 2 |V(^ 



(2) 
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where g is the SU(2)l coupling constant, rriq is the de- 
generate squark mass, and N X i is the projection of the 
lightest neutralino on the interaction eigenstates b, w°, 
hi, and h^- When the gaugino fraction approaches unity, 
we see that the decay rate is 



r = 
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(3) 



When the gaugino fraction approaches zero, there is an 
additional suppression from the light quark masses, 
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V. ANTIPROTON PRODUCTION 

Antiprotons are produced in the jets formed by the 
final state quarks. We are interested in the energy spec- 
trum of antiprotons dNp/dE at energies of a few GeV. 

It is quite difficult to calculate the spectrum dNp/dE, 
especially at low antiproton energies. There are basi- 
cally two theoretical frameworks available: a perturba- 
tive QCD approach and the Lund Monte-Carlo. In fig. |l| 



we compare the theoretical curves with the measured pro- 
ton spectra per quark jet in e + e~ annihilations at var- 
ious jet energies (14.5 GeV @, 29 GeV jll), and 45.6 
GeV jlH) and also show the theoretical predictions at 1 
TeV, one of the highest energies we need. 

In the perturbative QCD approach, Dokshitzer et al. 
Jl3| provide a formalism for calculating hadron spectra 
from parton spectra in the modified double-logarithmic 
approximation (MLLA). The full MLLA formula for in- 
clusive particle production is too difficult to evaluate, 
and no numerical expression is available. It has been 
shown that the simpler limiting spectrum reproduces the 
total charged particle spectrum quite well [Q. How- 
ever, as can be seen in fig. Ill neither the limiting spec- 
trum (long-dashed curves [|13|1 ) nor the low energy limit 
of the full MLLA spectrum (short-dashed curves Ej]) fit 
data for protons in a satisfactory manner. As a theoreti- 
cal alternative to the MLLA formalism, the Lund Monte 
Carlo jl5| has been used by several authors, most recently 
by Bottino ct al. p6[ . However, it also does not match 
the data for protons satisfactorily (histograms in fig. [j]). 

Inspecting the experimental data and the results of 
the Lund Monte Carlo, we observe that the (anti)proton 
production rate in the high-energy tail of the spectrum 
increases with jet energy. This is verified in the data 
for charged particle production in jets and is a generic 
prediction of the MLLA formalism Mjuj . 

Based on these arguments, we decide to take the value 
dNp/dE = 0.05 antiprotons /jet /GeV at the momentum 
p = 1.573 GeV in which we want the antiproton spectra 
(see sect. VII). This underestimates the production rate 
of antiprotons, and gives a conservative bound on the 
R-parity violating couplings. 

The number of jets per decay Nj Cts is equal to the num- 
ber of final state quarks, except when there is a top quark. 
The top decays before it can hadronize. We approximate 
the process t — > W + b using the fact that the low energy 
proton spectrum is insensitive to the jet energy. The bot- 
tom quark forms a jet. In addition, we take the hadronic 
branching fraction of the W , B(W — » hadrons) = 0.679, 
as the probability to form two additional jets in the de- 
cay W + — * ud and its charge conjugate. The final values 
for Vjets are 3 for a U C D C D C decay without top, 4.358 
for a U C D C D C with top, 2 for an LQD C without top and 
3.358 for an LQD C with top. 

We now obtain the volume production rate of an- 
tiprotons by multiplying the flat antiproton spectrum by 
the number density of neutralinos p x /m x and the num- 
ber of jets per decay Vj e t s , and dividing by the lifetime 

r x = i/r, 



dQp dNp p x A/j 



jets 



dE 



dE 



(5) 
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VI. ANTIPROTON PROPAGATION IN THE 
GALAXY 

We use the diffusion model of Webber et al. [|l7j as 
implemented by Chardonnet et al. (jlj| to describe dif- 
fusion of antiprotons through the galaxy. The flux of 
antiprotons at the outer solar system is approximated by 
multiplying the source spectrum by the antiproton ve- 
locity and an appropriate diffusion time. The diffusion 
constant in the galaxy depends on the rigidity p of the 
particle and can be approximated by the following ex- 
pression. 



K = 6 x 10 27 ( 1 



P 



3ZGV 



0.6 



(6) 



The region of turbulent magnetic fields responsible for 
diffusion can be approximated as a cylinder about 20kpc 
in radius and 2kpc high. This gives a diffusion time for 
antiprotons of 



t d = 6 x 10 15 (l 



P 



3GV 
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(7) 



The flux at the outer solar system beyond the influence 
of the solar wind is now given by 



4>t 



1 dQp 
4^~dE~ 



Vpt d - 



(8) 



This flux is in units of antiprotons s _1 cm -2 GeV -1 sr _1 . 

The low energy spectrum is affected by the solar wind. 
We use a simple model to account for this |l9[) . The 
observed antiproton energy E and momentum k are given 
in terms of the initial energy outside the solar system Ep. 



E7; = k r In 

P k c + 

E ¥ = E + AE, 



E 



E, 



AE, k < k c , 



(9) 



We adopt values of the critical momentum k c and en- 
ergy shift AE that correspond to the period of minimum 
activity in the 11 year solar cycle, k c =1.105 GeV and 
A_E=495 MeV. This gives the highest inner solar system 
flux, thus the conservative bound on the couplings. 
The solar modulated flux is now given by 



E 2 



El 

v 



(10) 



the top of the atmosphere is 6.4 x 10 _7 [cm 2 s sr GeV] -1 . 
We denote this flux as ^ ^ S (E). From the values of the 
observed energy E and momentum k, we determine the 
unmodulated p energy Ep. This is the energy at which 
we evaluate the diffusion time td and the velocity v. 

We now insist that the flux of antiprotons from neu- 
tralino decays be less than the observed flux. Including 
corrections due to solar modulation, the statement is 



Vpta p x jVjcts dN ¥ E 2 
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dE E 2 

p 



< $obs(£)- 



(11) 



We choose not to make a correction for the spallation 
production of antiprotons. The data are consistent with 
the picture that all cosmic ray antiprotons are made in 
spallation reactions, but there are a lot of uncertainties. 
There is no good model for the amount we should sub- 
tract, and we have chosen to ignore it. We thus obtain a 
conservative bound on the couplings. 

Using the data we have chosen, we obtain a model inde- 
pendent bound on the lifetime of a hadronically decaying 



relic with mass 



> 10 GeV 



> T p = 7.9 x 10 



28 , 



AW; 
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where we choose the mass fraction of neutralinos in the 
galactic dark halo as f x — Q, x h 2 . Since we assume that 
the neutralino is in the galactic halo at present, this 
bound applies for lifetimes longer than the age of the 
Universe i - I n other words, we can exclude lifetimes in 
the range t < t x < provided to < r£. 

For these long lifetimes, our bounds are in gen- 
eral better than those from diffuse photons J22], t x > 
7 x 10 25 s N- ]ets Vl x h 2 , and diffuse neutrinos p3fl , t x > 
10 25 s (m x /GeV) 17 / x , except for some supersymmetric 
models with neutralino masses in the TeV range. 

We can limit each element of the two coupling matrices 
A' and A" separately, by choosing the coupling such that 
all of the measured p flux comes from the neutralino de- 
cay via a single channel. Our excluded range of lifetimes 
then translates into the following excluded range of A 



Ap < |A| < A , 



where 



A ^ 3 - 48Xl °- 15 fe& 



(13) 



(14) 



VII. RESULTS 



is the upper limit on A coming from the antiproton flux 
and 



We proceed as follows. We choose a data point for the 
p flux at the earth [|fj. These data have large error bars, 
but will improve with the next generation of experiments 
to be done in space pl[ . According to the data from 
BESS, at an energy of 1.4 GeV the antiproton flux at 



A ~ 1.8 x 10~ 9 \/f (15) 

is the value of A for which the neutralino lifetime equals 
the age of the universe, which sets the validity of our 
analysis. In the previous formulas, m x is in GeV, f — 
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|A| 2 r x is in seconds, and the age of the universe is taken 
to be 10 10 yr. 

For virtually all models we find that Ap is smaller than 
Ao by at least three orders of magnitude. Conservatively 
taking -/Vj 0t s — 2 we find for the excluded range that 



£ = 7.3xl0 5 
Ap y m x /Gev 



x > 600. 



(16) 



If we would further demand that the neutralino be cos- 
mologically interesting, namely fl x h 2 > 0.025, we would 
expand the excluded range to Xo/Xp > 4000. 

We have plotted the upper bound on A'/ 12 as a function 
of the neutralino mass for 10 supersymmetric models in 
figure ||. The results are similar for all generations in 
the R-parity violating couplings, except for Ag jfc which 
shows the top quark threshold (see figure ||) . When the 
neutralino is lighter than the top quark, there is no bound 
on X'^ k at tree level. A bound can however be found at 
the one loop level Q . 

In our sample of supersymmetric models, we find the 
absolute upper limits 



< 2 x 10 



|A^- fc |p< 1 x 10 



-18 



-18 



ij3\P 
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< 3 x 10~ 19 



|A«*Ip<4x 10 



I A' 



I4 fc |p<lx 10 
|A" 3 Ip<1x 10 
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The apparent absence of an absolute upper bound on 
Ag'y is due to our neglect of one-loop diagrams. These 
upper limits depend on the neutralino mass, and improve 
considerably with higher masses. We find the following 
power law fits 



3.1 

l 

2 

2.7 
7 



(i,j,k^ 3) 
(j,k^3) 
(k ± 3) 

(i,k?3) 
(i ± 3). 



(23) 
(24) 
(25) 
(26) 
(27) 
(28) 



|A^ fe |p<6x 10- 22 (m x /TeV)- 
|A^ ife | P -<3x 10- 22 (m x /TeV)~ 3 
|Ak-lp<8x 10- 23 (m x /TeV)- 3 ' 
|A' ij3 lp<lx 10- 22 (m x /TeV)- 30 
|A^.| P -<3x 10- 22 (m x /TeV)- 2 ' 
|A-- 3 | P -<3x 10" 23 (m x /TeV)- 2 ' 

For comparison, the most stringent laboratory con- 
straints, those from proton decay p5j , bound not |A'| 
and | A" | separately but their product to |A'A"| < 
10 _24 (mo/TeV) 2 , where mo is the sfermion mass scale. 
These bounds are six orders of magnitude weaker than 
ours. On the other hand, because our analysis only ap- 
plies to lifetimes longer than the age of the Universe, our 
excluded region spans only 3 orders of magnitude. As 
mentioned in the introduction, the remaining gap of 3 
orders of magnitude may be filled by other cosmological 
considerations. 



We finally comment on the dependence of our bounds 
on model parameters. In figure ^ we plot the dependence 
on the squark mass scale m . As expected, the bound is 
better for lighter squarks. Regarding the dependence on 
neutralino composition, in general the bounds are bet- 
ter in the gaugino region. In the higgsino region, the 
bounds improve with heavier generations as the Yukawa 
couplings become larger. The dependence on other model 
parameters is negligible in comparison. 



VIII. CONCLUSIONS 

From observations of cosmic ray antiprotons we obtain 
strong constraints on each of the hadronic R-parity vio- 
lating couplings | A' | and |A"| that appear in extensions of 
the minimal supersymmetric model. Our analysis applies 
strictly in cases where the neutralino lifetime is longer 
than the age of the universe, and enables us to exclude 
a range Ap < |A| < Ao, where Ap is the upper limit from 
the antiproton flux and Ao is the validity boundary com- 
ing from the age of the universe. Using a large selec- 
tion of model parameters, we find absolute upper limits 

\Kjk\p ^ 2 x 10-18 and \K'jk\p ^ 10-18 when there is n0 

top quark in the final state, and |A^ 3fe |p < 4 x 10 -19 when 
there is. In virtually all models, we find that Ao > 10 3 Ap. 
Our bounds lie many orders of magnitude below the lab- 
oratory constraints coming from proton decay, and apply 
to each R-parity violating coupling separately. 
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FIG. 1. Antiproton spectra in quark jets from theoretical predictions and e + e annihilation data. Short- and long-dashed 
curves are the low-energy and the limiting MLLA spectra, histograms are obtained with the Lund Monte-Carlo. Data points 
are from TPC (14.5 GeV), TOPAZ (29 GeV), ALEPH (41.6 GeV, full squares), DELPHI (41.6 GeV, diamonds) and OPAL 
(41.6 GeV, open squares). The cross indicates our choice of dN^/dE at p — 1.573 GeV. 

FIG. 2. Bounds on A" 12 from cosmic ray antiprotons as a function of neutralino mass for 10 4 SUSY models. The scatter 
plots for \"j k are similar when i =fc 3. The scatter plots for A^ fc are also similar. 

FIG. 3. Bounds on A312 from cosmic ray antiprotons as a function of the neutralino mass for 10 4 SUSY models. The threshold 
effect is clearly seen. 

FIG. 4. Bounds on A'/ 12 from cosmic ray antiprotons as a function of the squark mass scale for 10 4 SUSY models. 
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Fig. g, E. A. Baltz and P. Gondolo, Phys. Rev. D 
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Fig. | E. A. Baltz and P. Gondolo, Phys. Rev. D 



O 



10 
10 
10 



-17 



-18 _ 



< 10 
10 



19 



■20 



-21 _ 



1 1 1 — I I I I I 



'"•f \-' i r.\ 



1 1 1 1 I I I I 



5 10 
I 10 

Oh 

p 

10 



-22 _ 



-23 _ 



-24 



10 



-25 



10 



-26 



10 



J l l 



J l l 



100 

m x (GeV) 



1000 







Fig. | E. A. Baltz and P. Gondolo, Phys. Rev. D 
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